      function qsat (t,p)
c          humidade especifica de saturacao
c
      implicit double precision (a-h,o-z)
      rasrv = 287.05/461.51
      etvq = 1. - rasrv
      espf=esat(t)/p
      qsat=rasrv*espf/(1.-etvq*espf)
      end


